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The visualization of chiral p-wave superfluidity in Fermi gases near p-wave Feshbach resonances 
is theoretically examined. It is proposed that the superfluidity becomes detectable in the entire 
BCS-BEC regimes through (i) vortex visualization by the density depletion inside the vortex core 
and (ii) intrinsic angular momentum in vortex free states. It is revealed that both (i) and (ii) are 
closely connected with the Majorana zero energy mode of the vortex core and the edge mode, which 
survive until the strong coupling BCS regime is approached from the weak coupling limit and vanish 
in the EEC regime. 
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There has been intense interest on p-wave pairing su- 
perfiuids recently in ultracold atom systems both exper- 
imentally and theoretically. It is highly expected that 
p-wave superfluidity may be realized experimentally, by 
using p-wave Feshbach resonance [1] which is now found 
for several atomic species ^Li and ^^K under magnetic 
field sweeping. In the two-dimensional (2D) system with 
the presence of magnetic field, the superfiuid has a spin- 
less chiral p-wave symmetry [H, Q . A theoretical task to 
be urgently clarified is how to detect the superfluidity. 
In the s-wave pairing case, the vortex imaging through 
the density depletion at a vortex core is decisive to es- 
tablish it [4]. There, the density depletion at vortex core 
is caused by the spontaneous particle-hole asymmetry Q 
of the discrete core- localized states, i.e^ so-called Caroli- 
de Gennes-Matricon (CdGM) states \&\. It is not quite 
obvious that the same is true for p-wave pairing 0, 0, , 
whose detailed study is one of the purposes here. 

The CdGM states in chiral p-wave superfiuid are in- 
triguing on its own right because it has exact zero energy 
state (ZES) as shown in the BCS limit where the chemi- 
cal potential fi is equal to the Fermi energy Ep 
Since the ZES can be described by the one-dimensional 
(ID) Majorana equation [10], the quasiparticle of the 
ZES is a Majorana fermion. Different from ZES of the 
Jackiw-Rebbi solution in the ID Dirac equation at the 
domain wall [jjj, ll2|, ll3| , the remarkable feature of Ma- 
jorana fermions is that their creation is expressed by the 
self-Hermitian operator, and their host vortices obey the 
non-abelian statistic s [1 , These can be utilized for 
quantum computer [15|. It is natural to expect that 
the observation of p-wave Feshbach resonance [l| opens 
a door connecting ultracold atoms to diverse research 
fields, such as quantum computation, the condensed mat- 
ter, and the relativistic field theory. 

The other low-lying mode characteristic to p-wave su- 
perfiuid is the edge mode, resulting from the Andreev 
scattering at the rigid wall of the system Due to the 
edge mode, the p-wave superfiuid has the spontaneous 
mass flow in the vortex free state, which accompanies 



the angular momentum arising from the internal motion 
of p-wave pairs [17]. The observation of macroscopic an- 
gular momentum could be used as another method to 
detect evidence of the p-wave superfluidity. 

The merit of study in atomic gases is that we can 
manipulate the interaction through Feshbach resonance. 
The Cooper pairs realized within /i > are turned to 
strong coupling at /i ^ +0, and further to Bose-Einstein 
condensation (BEG) of molecular bosons when /i < 0. 
Contrasted to the BCS-BEC crossover in the s-wave pair- 
ing, there is a topological phase transition at /i = be- 
tween BCS and BEG in the p-wave pairing [14]. The 
purpose of this Letter is to examine the possibility to 
experimentally detect the evidence of superfluidity as- 
sociated with p-wave Feshbach resonances, in the entire 
range from BCS to BEG regimes. It is proposed that the 
superfluidity becomes detectable through the vortex vi- 
sualization and intrinsic angular momentum. From fully 
microscopic point of view based on the Bogoliubov-de 
Gennes (BdG) equation, which enables the systematic 
study in the BCS-BEC regimes, we reveal that their vi- 
sualization is closely connected with novel quasiparticle 
states, such as the Majorana ZES and the edge mode. 

The attractive p-wave interaction by the Feshbach res- 
onance can produce a pairing between spinless fermions 
[1, 0] . The pair potential for chiral p-wave pairing in 2D 



systems can be expressed as A(r, /c) = 



A+i(r). 



-A_i(r), in terms of the eigenstate of the orbital 



angular momentum Lz\m) = m\m) (m = 0,±1), where 
r and k are the center-of-mass coordinate and the rel- 
ative wave vector of the pair. Under the pair potential 
A(r, k)^ the quasiparticle eigenstates with the wave func- 
tion [uq^VqY is described by the BdG equation 0, [l^]. 
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where h = ks = l- The applicability of the sin- 
gle channel model in the p-wave case is discussed in 
Ref. [l^, where it is verified that this simple theory is 
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suitable for describing the BCS-to-BEC regimes associ- 
ated with narrow p-wave Feshbach resonances. We set 
n(r) = -^E±{A±i(r)n± + ^n±A±i(r)}andi7oW = 

+ V{r) - ji with mass M, kp = \/2MEf, and 
the harmonic trap potential V{r) which confines atomic 
gases. While we mainly discuss the case when V{r) =0, 
we also confirmed that the presence of a harmonic trap 
does not change essential features obtained with V = 0. In 
the 2D polar coordinates r = (r, [I\± = ^^^{dr^^de) . 
Assuming the cylindrical symmetry of the pair potential, 
Am(^) = A^(r)e*'^'^^ with the winding number Wm ^ ^• 
Thus, the eigenstates labeled by q in Eq. ([1]) are char- 
acterized by the azimuthal quantum number g^i G Z as 
Uq(r) = Uq{r)e'^^^, Vq(r) = Vq(r)e^^^^-^+i-^)^, and we 
find W-i =W-^i + 2. In this work, we consider the case 
when the dominant chiral component is A+i(r) and has 
vortex winding W-^i = — 1, i.e., vortex current fiows along 
the direction opposite to the orbital motion of the chiral 
pairing. The passive component A_i(r) are only induced 
around the vortex core and the boundary. The other 
cases of vortex winding will be discussed in Ref. [18]. 

The pair potentials A±i are determined by the 
gap equation A±i(r) = ^ E£;^<o[^q(^)°T^q(^) - 
Uq{r)\3z^Vq{r)] at zero temperature (T = 0) [9]. From 
practical point of view, it is convenient to express the 
coupling constant g with the two-body bound state en- 
ergy in vacuum E,, as g'^ = - j j^^^f^^ g 0. 
Note that in contrast to the 5- wave case, the resulting 
gap equation still contains the logarithmic divergence on 
the energy cutoff Ec [20]. Throughout this work, we cal- 
culate the set of equations with Ec = 20, 40, and QOEp 
and verify that there is no qualitative difference between 
them. Hence, in this Letter, we present the results with 
Ec = 60Ef- We solve the set of equations in a cylinder 
with the radius R = 50kp^ and the boundary condition 
Uq{r = R) = Vq{r = R) = 0. Here, all the quantities 
are scaled by the energy unit Ejr and the length unit 
kp^. During the self-consistent calculation, we vary /i so 
as to fix the total particle number N = J p{r)dr = 604, 
with the particle density p{r) = <o When 
changing the pairing interaction by varying £^5 within 
— 1.8<^5/£^F<l-0, the chemical potential and the maxi- 
mum of the pair potential A = max |A+i(r)| continuously 
change from ii = {).^Ef and A = 0.79£;f at Ei, = 1.{)Ef of 
strong coupling BCS to ii = -{).21Ef and A = 1.46£^f at 
Ei, = -I.^Ef of EEC. 

In Fig. [H^a), we plot the energy spectrum of the vor- 
tex state at £^5 = 1.0£^f- There is a symmetry Eq^ = 
—E-q^. Two branches appear inside the bulk energy 
gap ^gap = O.AEf' The branch labeled as "Edge" in 
Fig. [Tfa) consists of the eigenstates with a wave function 
localized at the edge of the cylinder as the Andreev bound 
state. The other branch labeled as "CdGM" is composed 
by the core-localized mode. The energy separations be- 
tween discrete levels of the CdGM states are large in this 
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FIG. 1: (Color online) (a) Energy spectrum Eq of the = 
— 1 vortex as functions of qe^ at Eh — I.^Ef with fi^O.GEp 
and A = 0.79^F. The dashed and dashed-dotted lines denote 
Eq^"^ and Eq^'^ . (b) Low- lying eigenenergies at qe = in the 
w-\-i = — 1 vortex as functions of Et. Circles and triangles 
denote the eigenenergies of the ZES and the lowest bulk ex- 
citation state, respectively. The dotted-dashed line denotes 



strong coupling. In the BCS limit, dispersion relations of 
the branches can be analytically derived from Eq. ^ as 

E^"^ = —{qe — ^^2^"*" ) for edge mode by Stone and 

Roy [16|, and = -{qe - ^^)^o + (n - 
withnGZ, w^i^O and cjq ~ A^/^;^ <Ccji A for CdGM 
mode using \qQ\ <C /c^?^ and A <C /i ~ Ef by us [ist . 
These analytical expressions in the BCS limit, presented 
by lines in Fig. [ija), still give nice fitting to our numer- 
ical results of the strong coupling. It is noted that the 
eigenenergies of E^"^ and E^''^ can become exactly zero 
in vortices with arbitrary odd number of W-^i. 

We notice that the behaviors of the -dependence of 
ZES and the density depletion can be classified to three 
regimes: (i) the BCS limit of weak coupling with (ic^Ef 
(El) I Ef ^ 1) , (ii) the resonance regime of strong coupling 
BCS with < /i<^Ef {\Eb\/EF < 1.0), (iii) the BEG 
regime with /i<0 {Ei)<—1.0Ef)- The topological phase 
transition occurs at /i = (£^5 = —I.OEf)- 

The low- lying eigenenergies of the spectrum with qe = 
are displayed in Fig.llfb) as a function of £5. The lowest 
extended states denoted by triangle points correspond to 
the bulk excitation gap £gap- In the resonance regime 
(ii), £gap ^ M £gap approaching the topological 
phase transition in the p-wave pairing. The bulk excita- 
tion in the BEG regime (iii) is again gapful as £gap ~ — /i, 
corresponding to the binding energy of molecules. Within 
the excitation gap, doubly degenerate ZES's of CdGM 
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and edge modes behave as follows. In the BCS limit, it 



is known from E^"^ and the index theorem 13, 21 1 that 



the ZES always appears. Figure [Tfb) shows that the ZES 
survives also in the resonance regime beyond the BCS 
limit, until just before the topological phase transition 
at Ei) = —1.0Ef. We have also confirmed the degenerate 
ZES's appear in harmonically trapped systems [l8|. In 
the EEC regime beyond the topological transition, the 
ZES vanishes and merges to the excitation gap. In the 
vicinity of £^6 = — 1.0£^f, the ZES shifts to finite energies 
by the interference of the CdGM and the edge states. The 
interference of the ZES's becomes effective when /i +0, 
because the length scale of the spatial variation of Uq and 
Vq becomes comparable with the system size R [22]. This 
shift of ZES is seen at Ei, = —I.OEf in Fig. [T^b), since 
this change occurs slightly before Ei) = — I.OEf when the 
system radius R is finite. 

The quantized vortex which is a hallmark of superflu- 
idity is detected by the depletion of the particle density 
p{r) at the vortex core. Here, we discuss the relation 
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FIG. 2: (Color online) (a) Particle density p(r) in the w-^i = 
— 1 vortex state at Eb = I.OEf (solid line) and the BEC regime 
Eb = —I.SEf (dotted-dashed line). The dashed line shows the 
contribution of the extended states at \Eq \ >0.4:Ef to p{r) at 
Eb^l.OEF, and the shaded area depicts the contribution of 
the ZES. The dotted line denotes constant p(r) in the BCS 
limit The inset in (a) shows the density depletion at 

the core p{r = 0)/p(oo) (triangles) and the contribution of 
the ZES (circles), (b) Local density of states N'{r = 0,E) at 
Eb/Ep = 1.0, -0.4, and -1.8. 



of the ZES and vortex visualization by p{r = 0) at the 
core. To see the contribution of the ZES and the ex- 
tended states outside the gap, in Fig. Efb) we present 
the local density of states (LDOS) at r = defined as 
^/{r,E) = \uq{r)\^5{E - Eq). In the BCS limit, it 
is known that p(r) keeps constant even inside the vortex 
core when Wj^i = — 1 [9], as presented by dotted line in 
Fig.[2l[a). This indicates that the vortex is completely in- 
visible through the density imaging in the BCS limit. As 
approaching the resonance regime from the BCS limit, 
the vortex core gradually becomes visible by the contri- 
bution of extended states. We show p(r) at £^5 = I.OEf 
by solid line in Fig. [2fa). If the contributions of CdGM 
states are removed from p(r), the density depletion by 
the extended states is much deeper as shown by dashed 
line in Fig. EJ^a), which reflects the spatial structure of 
A+i(r). This strong depletion at the vortex core is com- 
pensated by the ZES contribution, as shown by the shade 
in Fig.[2ja), and it makes the depletion weaker. This con- 
tribution of ZES is absent in the case of 5-wave pairing. 
Thus, the density contrast of vortex core visualization 
is weaker in the chiral p-wave pairing. The vortex vis- 
ibility p{r = 0)/p{r = 00) and the contribution of the 
ZES are plotted in the inset of Fig. [21(a) as a function 
of Eb. When — 1 < E^^/Ef < in the resonance regime 
(ii), p(r = 0) comes from only the ZES, since the LDOS 
at Ei) = —O.AEf in Fig. [2](b) has no intensity outside the 
gap for Eq<0. Note that the absence of the negative en- 
ergy band results from the strong coupling effects, such as 
\p\ <^ A ^ Ef, while there still remains the particle- hole 
symmetry. In the BEC regime, the vortex core is clearly 
visible through the density imaging experiments [3]. As 
shown in Fig. [2](b) for £^5 = — 1.8£^f, the ZES disappears 
and the excitation spectrum becomes gapful. Thus, as 
the dotted-dashed line in Fig. [21(a), p{r) approaches per- 
fect depletion p(r) ex |A+i(r)p in the BEC. 

Next, we discuss the relation of energy spectrum and 
intrinsic angular momentum in the vortex free state. In 
Fig. [31(a) we plot the energy spectra of the vortex free 
state with W-^i = ai E^ = —O.SEf in the resonance 
regime. There we see the linear branch of the edge mode 
at low energies within the gap. In the vortex free state, 
the edge mode has small but finite energy at ^51 = 0, dif- 
ferent from the vortex states. The edge states with neg- 
ative energies are occupied at T = 0. Due to asymmetry 
on qe, the occupied edge modes with qe > produce the 
spontaneous mass current [16|] je{r) = ^\uq{r)\'^ 
around the edge of cylinder. Hence, net angular mo- 
mentum Lz = J rjo{r)dr is not zero even in the vortex 
free state. We numerically find that L^/N = 0.509 at 
Ei) = — O.SEf. This is also understandable as the angular 
momentum carried by the orbital motion of the Cooper 
pair with the chirality m = +l. The small deviation from 
the intuitive expectation +^ is due to the pairing A_i(r) 
with W-i =+2 induced around the edge. 

More noteworthy is that there still exists the spon- 
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FIG. 3: (Color online) Energy spectrum Eq of the vortex free 
state near the resonance: (a) Eb = —O.SEf with /i = 0.07^F 
and A/Ef = 1.29 and (b) Eb = -1.2Ef with fi = -0.05Ef 
and A/Ef = 1.36. The inset in (b) depicts the current density 
\je{r)\ at Eb = -l.2EF within |x|, |^| < 50A;^\ 

taneous mass flow even in the BEC regime with /i < 0. 
Figure mb) shows that the branch of the edge mode with 
the positive (negative) eigenenergy merges into the upper 
(lower) band of the bulk excitations, and the resulting en- 
ergy spectrum becomes gapful with the gap Hence, 
the whole spectrum remains asymmetry with respect to 
qo^ which can produce the spontaneous mass flow along 
the edge, even without edge modes in the BEC regime. 
The inset in Fig. [Sjb) shows je{r)^ which is localized at 
the edge within It is also found that Lz/N stays 

around 0.5 within 15% in the range — l.S< Ei)/ Ejr < 1.0. 

Finally, we examine influences of a trap potential V{r) 
on the edge mode and the mass flow. Figure (H^a) de- 
picts the energy spectrum of N = 1,000 atoms in the 
vortex free state under a harmonic trap V{r) = ^mcj^r^, 
which still accompanies low- lying edge modes. The pres- 
ence of the edge mode in a harmonic trap was dis- 
cussed in Ref. [23[ under a given uniform pairing field 
A+i(r) = const. We examine the edge mode by the full 
self-consistent calculation including spatial variation of 
A±i(r). In the presence of F(r), since the wave functions 
of the edge mode have a broad peak around r~0.85i?TF 
with Thomas- Fermi radius i^xF, the spontaneous cur- 
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FIG. 4: (Color online) (a) Energy spectrum Eq of the vortex 
free state with w-\-i =0 in harmonically trapped state, where 
fi = 0.78EF and max | A+i (r)| = 1.05^f with ^f/cj = 42.5 at 
^^ = 0.0. The solid line denotes E^q""^ with i? = i?TF = 74.5A:^^ 
and A=|A+i(r = 0.85i?TF)|=0.35^F. (b) The corresponding 
current density |i6'(^)| within |x|, |y| < 1.05i?TF. 



rent je{r) has broad profile around Rtf^ as shown in 
Fig. Sl^b). However, this does not seriously change the 
result 

The azimuthal flow along the Thomas-Fermi edge 
may affect the collective surface excitations, such as the 
quadrupole mode, which gives rise to the frequency shift 
between the co- and counter-rotating surface motion [IJ] • 
Hence, the measurement of the intrinsic angular momen- 
tum Lz/N 0.5, through the shift of the quadrupole 
frequencies, provides evidence of p-wave superfluidity. 

In summary, we have theoretically investigated novel 
features of chiral superfluid Fermi gases in the strong 
coupling BCS regime and the BEC regime near p-wave 
Feshbach resonance. Majorana zero energy states of the 
vortex core mode and the edge modes survives from the 
BCS limit {fi ^ Ep) to the BCS-BEC transition point 
(/i = 0), and vanish in the BEC regime. In the relation to 
the low-energy modes, we estimate (i) vortex visualiza- 
tion by the density depletion and (ii) the intrinsic angular 
momentum in the vortex free state, finding that these can 
be used to experimentally detect the appearance of su- 
perfluidity in wide range from BCS to BEC through the 
p-wave Feshbach resonance. 
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